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We study the general features of the dynamics of the phantom field in the cosmological context. In 
the case of inverse coshyperbolic potential, we demonstrate that the phantom field can successfully 
drive the observed current accelerated expansion of the universe with the equation of state parameter 
w<f, < — 1. The de-Sitter universe turns out to be the late time attractor of the model. The main 
features of the dynamics are independent of the initial conditions and the parameters of the model. 
The model fits the supernova data very well, allowing for —2.4 < < —1 at 95 % confidence level. 



PACS numbers: 98.80.Cq, 98.80.Hw, 04.50.+h 
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I. INTRODUCTION 

The observational evidence for accelerating universe 
has been for the past few years one of the central 
themes of modern cosmology. The explanation of these 
observations in the framework of the standard cosmol- 
ogy requires an exotic form of energy which violates 
the strong energy condition. A variety of scalar field 
models have been conjectured for this purpose including 
quintessence El, K-essence Q and recently tachyonic 
scalar neldsHH 1 1 Q . A different approach to the 
cosmic acceleration is advocated in Rcfs |3.(For a review 
of the issues related to the cosmological constant and 
the dark energy, see EHH El 13). By choosing an 
appropriate potential V(4>) and tuning the parameters of 
the model one can account for the current acceleration 
of the universe with Q$ = 0.7 and f2 OT = 0.3. However, 
all these models lead to the equation of state parameter 
> — 1. The recent observations seem to favor the 
values of this parameter less than -1. 

A scalar field with negative kinetic energy called the 
phantom field is proposed by Caldwell, Carroll et al and 
others to real ize the possibility oblate ^ me accele ration 

a field has a very unusual dynamics as it violates null 
dominant energy condition (NDEC). The models with 
equation of state parameter less than -1 are known to 
face the problem of future curvature singularity which 
can however be overcome in specific models of phan- 
tom field. Inspite of the fact that the field theory of 
phantom fields does encounter the problem of stability 
which one could try to bypass by assuming them to be 
effective fields pa . l3lt l33j | . it is nevertheless interesting to 
study their cosmological implications. Curiously, phan- 
tom fields have glorious lineage in Hoyle's version of the 



Steady State Theory. In adherence to the Perfect Cos- 
mological Principle, a creation field (C-field) was for the 
first time introduced [11] to reconcile with homogeneous 
density by creation of new matter in the voids caused by 
the expansion of the universe. It was further refined and 
reformulated in the Hoyle and Narlikar theory of gravita- 
tion [H, see RefsHlUllillil f ° r details ° f the C " field 
cosmology. The C-field appeared on the right hand side 

of the Einstein equation. What was conserved was the 
sum of the stress-tensors of matter and C-field and nei- 
ther was conserved separately. The C-field violated the 
weak energy condition. In this context, it may also be 
recalled that wormholes do require violation of energy 
conditions, in particular of the average null energy con- 
dition m| . Interestingly, it has recently been shown that 
a very finite amount of exotic energy condition violating 
matter is sufficient to produce a wormhole |4o| . Thus 
phantom fields though very exotic are not entirely new 
and out of place. Here the motivation is strongly driven 
by the observation. 



II. DYNAMICS OF THE PHANTOM FIELD 

The Lagrangian of the phantom field minimally cou- 
pled to gravity and matter sources is given by[31| 



£= (WnG)- 1 R+-g^d^ 



,<t>-V(<t>)+C s 



(1) 



where C source is the remaining source term (mat- 
ter, radiation) and V(4>) is the phantom potential. The 
kinetic energy term of the phantom field in Q enters 
with the opposite sign in contrast to the ordinary scalar 
field (we employ the metric signature, -,+,+,+). It is the 
negative kinetic energy that distinguishes phantom fields 
from the ordinary fields. The Einstein equations which 
follow from (JTJ are 
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where Tj£ and T^° urce are the stress tensors of phan- 
tom field and of the background (matter, radiation). In 
contrast to the case of the creation field for which the 
individual stress-tensors could not be conserved indepen- 
dently, the phantom field energy momentum tensor 
here and the usual matter field stress-tensor T^° urce are 
indeed conserved separately. The former was introduced 
to create new matter out of nothing so as to keep the den- 
sity uniform in an expanding universe. The latter has no 
such purpose instead its job is to accelerate expansion at 
late times with < — 1. 

The stress tensor for the phantom field which follows 
from (JTJ has the form 



rpph 



-du<f>du(, 



9n» 



(4) 



We shall assume that the phantom field is evolving in 
an isotropic and homogeneous space-time and that 4> is 
a function of time alone. The energy density p^ and 
pressure 



obtained from are 



+ V(<j>), P4> 



V(4>). (5) 



The fact that the phantom field feels, in contrast to ordi- 
nary field/matter, opposite curvature of the space-time 
is reflected in the negativity of kinetic energy in eqn. 
The equation of state parameter is now given by 



(6) 



Now the conditions w^, < 

A2 



1. 



, > would prescribe 
the range, < (p 2 < 2 V (</>). The Friedman equation 
which follows from eqn. J3J) with the modified energy 
momentum tensor given by eqn. J3J is 



a 2 {t) _ 8ttG 



[p<p + Pb] 



(7) 



where the background energy density due to matter and 
radiation is given by 



„ PR , Pm 
Pb = — + — T 

a 4 er 



(8) 



The peculiar nature of the phantom field also gets re- 
flected in its evolution equation which follows from eqn. 



4> + —4> = v,M- 

a 



(9) 



Note that the evolution equation © for the phantom 
field is same as that of the normal scalar field with the 
inverted potential allowing the field with zero initial ki- 
netic energy to roll up the hill; i.e. from lower value of 




FIG. 1: Evolution of the phantom field is shown for the model 
described by eqn. 1131 . Due to the unusual behavior, the 
phantom field, released with zero kinetic energy away from 
the origin, moves towards the top of the potential. It sets into 
the damped oscillations about = and ultimately settles 
there permanently. 



the potential to higher one. At the first look, such a sit- 
uation seems to be pathological. However, at present, 
the situation in cosmology is remarkably tolerant to any 
pathology if it can lead to a viable cosmological model. 

The conservation equation formally equivalent to eqn. 
@, has the usual form 

P4,+iH(jt^+P4,) = Q (10) 
and the evolution of energy density is given by 

P4> = J " (11) 

with 

1 



C(o) = 



(K e /P e 



1 



where the ratio of kinetic to potential energy (K e /P e ) is 
given by 



Kr 



(12) 



The evolution of potential to kinetic energy ratio plays 
a significant role in the growth or decay of the energy 
density at a given epoch and will be crucial in the 
following discussion |47j . 

Let us next address the question of choice of the poten- 
tial of the scalar field which would lead to a viable cos- 
mological model with w$ < — 1. An obvious restriction 
on the evolution is that the scalar field should survive till 
today (to account for the observed late time accelerated 
expansion) without interfering with the nucleosynthesis 
of the standard model and on the other hand should also 
avoid the future collapse of the universe. It was indicated 
by Carroll, Hoffman and Trodden [2(| as how to build 
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FIG. 2: Phase portrait (plot of Y 2 = 4>/Mp versus Yi = 
(f>/M p ) of the model described by eqn. I13H . Trajectories 
starting anywhere in the phase space end up at the stable 
critical point (0,0). 



models free from future singularity with < — 1. In 
this paper we examine a model with the scalar phantom 
field which leads to a viable cosmology with < — 1 
and shortly after driving the current acceleration of uni- 
verse the field settles at w$ = — 1 thereby avoiding the 
future singularity. We confront the model with super- 
nova la observations to constrain its parameters. The 
model fits the supernova data quite well for a large range 
of parameters. 



III. A VIABLE MODEL WITH THE EQUATION 

OF STATE, w<f, < -1 



We shall here consider a model with 



V(<j>) = V Q 



cosh 



a<j) 



~] -l 



(13) 



Due to its peculiar properties, the phantom field, released 
at a distance from the origin with zero kinetic energy, 
moves towards the top of the potential and crosses over 
to the other side and turns back to execute the damped 
oscillation about the maximum of the potential (see fig- 
ure After a certain period of time the motion ceases 
and the field settles on the top of the potential perma- 
nently to mimic the de-Sitter like behavior (w^ = —1). 
Indeed, the de-Sitter like phase is a late time attrac- 
tor of the model (see Fig. It should be emphasized 
that these are the general features of phantom dynam- 
ics which are valid for any bell shaped potential, say, 
V{4>) = Vq [l + a0 2 /Afp] or the Gaussian potential 
considered in Ref. [2q . 

In order to investigate the dynamics described by eqns. 
and |(3J), it would be convenient to cast these equations 
as a system of first order equations 



1 



H(Y 1 ,Y 2 ) 



Yn 



(14) 



FIG. 3: The ratio of kinetic to potential energy of the phan- 
tom field is plotted for the potential given by eqn. 1131 1 fol- 
ia — 2 (solid line) and a = 3 (dashed line). The evolution 
of \K e /P e \ starts later but peaks higher for larger value of a. 
The height of the peak is independent of Vq. The change in 
the value of Vo merely shifts the position of the peak. 



Y' 



SY 2 



H(Y U Y 2 ) 



dV{Y{ 



dYi 



where 



Yn 



M 2 '' 

p 



V 



V(Yi) 



(15) 
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and prime denotes the derivative with respect to the vari- 
able N = ln(a). The function H(Yi, Y 2 ) is given: 



H{Y U Y 2 



Y 2 



■V(Yi) 



Pb 
M 4 



(17) 



where p b = p l r e^ N + pi n e~ 3N . Using eqns. ©, d 
and lfl"5|) . it is not difficult to see that {Yi,Y 2 ) = (0,0) is 
a fixed point of the system. Numerical analysis confirms 
the stability of the fixed point (Fig. |2J). Thus the de- 
Sitter like solution is the late time attractor of the model. 

As for the initial conditions, for convenience we shall 
set them in the radiation dominated era with a; = 1 and 
pi = 1 MeV 4 . Note that at the present epoch, the scale 
factor a would nearly be 4 x 10 9 . The initial value for 
(j) as well as the values of parameters in the potential 
(a, Vo) are chosen so as to ensure a viable cosmological 
model. These choices could be construed as fine tuning 
in the model. But it is no worse than the fine tuning in 
some models of quintessence, for instance, see Ref . |47j. It 
should however be noted that quintessence models based 
on tracker potentials are independent of initial conditions 
and only require fine tuning of parameters of the poten- 
tial. 

We shall now describe the dynamics of the model due 
to the potential l|13fl . Initially, the field is displaced from 
the maximum of the potential and its energy density is 
subdominant with respect to the background energy den- 
sity. As a result the background plays the deciding role 
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FIG. 4: The energy density is plotted against the scale fac- 
tor: solid line corresponds to for a = 1.26 in case of the 
model CHI with V 1/4 ~ 3 x 10" 30 M P . The dashed and dot- 
ted lines correspond to energy density of radiation and mat- 
ter. Initially, the energy density of the phantom field is ex- 
tremely subdominant and remains to be so for most of the 
period of evolution. At late times, the field energy density 
catches up with the background, overtakes it and starts grow- 
ing (w^j < — 1) and drives the current accelerated expansion 
of the universe before freezing to a constant value equal to — 1 
(in future). 
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FIG. 6: Dimensionless density parameter SI is plotted against 
the scale factor for the model described by eqn. 1131 with 
a = 1.26 and V 1/4 ~ 3 x 10~ 30 Af p for: (i) phantom field 
(solid line), (ii) radiation (dotted line) and matter (dashed 
line). Late time behavior of the phantom field leads to the 
present day value of Q^, — 0.7 and Q m = 0.3 for the equation 
of state parameter w^, < — 1. 
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FIG. 5: Evolution of the equation of state parameter w^, is 
shown as a function of the scale factor for the model described 
in figure 0] in case of a — 2.5 (solid line) and a — 2 (dashed 
line) . Except for a short period, is seen to be constant 
(-1). The parameter evolves to negative values less than -1 
(smaller for larger value of a) at late times leading to cur- 
rent acceleration of the universe . It then executes damped 
oscillations and fast stabilizes to wj, = — 1. 



in the evolution dynamics. The Hubble damping due to 
Pb >> p<fi is (extremely) large in the field evolution equa- 
tion. Consequently, the field does not evolve and freezes 
at the initial position mimicking the cosmological con- 
stant like behavior. Meanwhile, the background energy 
density red shifts as 1/a™ (n=4 for radiation). The phan- 
tom field <f) continues in the state with w$ — — 1 till the 
moment p^ approaches pb- The background ceases now 
to play the leading role (becomes subdominant) and the 



phantom field takes over and it moves fast towards the 
top of the potential. The kinetic to potential energy ra- 
tio \K e /P e \ rapidly increases and becomes maximum at 
4> = allowing the equation of state parameter to at- 
tain the minimum (negative) value. This leads to the 
fast increase in p^. Damped oscillations then set in the 
system making the ratio to oscillate to zero (see Fig. |3J) 
and allowing p^ to settle ultimately at a constant value 
(w0 = —1) for ever. This development is summarized 
in Figs. I0J and (JSJ) . As shown in Fig. @, the ratio 
| .Kg /P e | takes off later but peaks higher for larger value 
of the parameter a. The reason is that for larger value 
of a the field energy density is lower and consequently 
it takes longer for p^ to catch up with the background. 
Once the background value is reached, the field sets into 
motion and rolls towards the maximum of the potential 
and the roll-up is faster as the potential gets steeper;i.e. 
larger is the value of a. We should however emphasize 
that the main features of the evolution are absolutely in- 
dependent of the initial conditions and the values of the 
parameters in the model. However, tuning of the param- 
eters is required to get the right things to happen at right 
time. For the major period of time the equation of state, 
w$ = — 1 while for relatively short time w$ is required to 
be < — 1 (see Fig. ©). The later happens when p$ over- 
takes the background (at late times) and starts growing, 
leading to the fast growth of f2^. By tuning the parame- 
ters of the model it is possible to account for the current 
accelerated expansion of the universe with 17^, = 0.7 and 
fl m — 0.3 during the period when < —1 ( see Fig. [JJ. 
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IV. CONSTRAINTS ON PARAMETER SPACE 
FROM SUPERNOVA OBSERVATIONS 

We use the Supernova la observations to put constrains 
on the parameter space of the phantom field model. For 
that let us first note that the luminosity distance (d/,) 
for a source at redshift z located at radial coordinate dis- 
tance r is given by = (1 + z) r, where clq is the 
present value of the scale factor. This can be used to de- 
fine the dimcnsionlcss luminosity distance, Pl = HqcIl, 
where Hq is the present value of the Hubble parameter. 
The apparent magnitude (m) of the source is given as 

m(z) = M + 5 log[P L (z)] (18) 

where M. = M — 5 \ogHo+ constant. 

We use the magnitude-redshift data of 57 supernovae 
of type la, which include 54 supernovae considered by 
Perlmutter et al (excluding 6 outliers from the full sample 
of 60 supernovae) [H, SN 1997ff at z = 1.755 [H and 
two newly discovered supernovae SN2002dc at z = 0.475 
and SN2002dd at z = 0.95 The best fitting values of 
the parameters can be obtained through x 2 minimization 
where, 

X ^ Smf ' 1 ' 

i—l L 

Here mf refers to the effective magnitude of the ith 
supernovae which has been corrected by the lightcurve 
width-luminosity correction, galactic extinction and the 
K-correction. The uncertainty in mf is denoted by 
Smf. 
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Redshift 

FIG. 7: The magnitude-redshift plot for 57 supernovae, in- 
cluding the recently discovered, SN2002dc and SN2002dd. 
The theoretical curve is for the best fit value obtained at 
a = 2.3 and fiy = 3V a /p c = 1500 for fa = 3.5Af P . 

Since our aim is to constrain parameter space of the 
phantom field we tuned the p m at initial epoch in such 



a way that in the absence of phantom dynamics it yields 
the ratio of matter density to critical density (p c ) today 
as fi m = 0.3 which is consistent with other observations 
including WMAP [H]]]. However, since the equations for 
the matter and radiation density evolution and the phan- 
tom dynamics are coupled, not all values of the parame- 
ters Vo and a would yield Q m = 0.3 and fi<£ = 0.7 today. 
Apart from these two parameters, there are initial condi- 
tions on <fi and <f> of which the latter was fixed to zero (in 
fact we found that results are not affected on reasonable 
variation of this initial condition). For a viable evolution 
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FIG. 8: The 68.3 % and 99.7 % confidence regions of the 
phantom parameter space are shown as black and grey regions 
respectively 

of the universe it is required that initial value of 4> 
be at least of the order of few. In case of <\>i much larger 
than unity, the equation of state turns out to be close to 
-1 for admissible values of fl m and Sl^, thus reasonable 
values of </>,; are of the order of few which yield present 
value of the equation of state parameter less than -1. The 
best fitting parameters after \ 2 minimization were found 
to be a = 2.3, fl Vo = 1500 {n Vo = 3V /p c ), <f) t = 3.5M P 
and M = 23.8 with \ 2 — 61.9 and x 2 P er degrees of 
freedom as 1.17 which represents a good fit, as is shown 
in Fig. 7. For the same data set the % 2 per degree of 
freedom for the best fitting flat model for the constant 
equation of state turns out to be 1.08 with O m = 0.32 
[52j | . The best fit phantom parameters yield an accelerat- 
ing universe with present values of the equation of state 
as w<p = —1.74, n m = 0.3 and fl<p = 0.7 

In order to constrain the parameter space we fixed the 
arbitrariness in (pi to its best fit value and marginalized 
over M.. In Fig. 8 we have shown the 68.3 % and 99.7 
% confidence regions in the parameter space of f2y and 
a. As depicted in the figure, a large region of the pa- 
rameter space is allowed by the supernova observations. 
However, we should emphasize that a small change in 
a corresponds to a large variation in f2vb which reflects 
a fine tuning of parameters similar to other models of 
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dark energy. Fig. 9 depicts 95.4 % confidence region of 
the parameter space with different allowed ranges of the 
equation of state. Dark energy models with the equation 
of state less than -1 have recently been analyzed and the 
bounds have been obtained il |H |27|,[35|. At 95.4 
% confidence level the bound on for the model un- 
der consideration is found to be -2.4 < w<p < -1. Thus, 
the phantom field model being constrained by supernova 
observations, favors a lower value of equation of state 
parameter. Similar bounds have been obtained for other 
models with phantom energy [24| . 




1 1.25 1.5 1.75 2 2.25 



a 

FIG. 9: The figure shows 95.4 % confidence region of the 
phantom parameter space. The light gray shaded region 
shows allowed parameters which yield > —1.3, the dark 
shaded gray region corresponds to —1.3 < w^, < —1.6, 
whereas the black region represents w$ < —1.6. 



V. DISCUSSION 

In this paper we have investigated the general features 
of the cosmological dynamics of the phantom field. In 
the case of the inverse coshyperbolic potential, we have 
shown that the phantom field can account for the cur- 
rent acceleration of universe with the negative values of 
the equation of state parameter {w^ < — 1). The general 
features of the model are shown to be independent of the 
initial conditions and the values of the parameters in the 
model. However, tuning of the parameters is necessary 
for a viable evolution. Unlike, the quintessence mod- 
els based on tracker potentials which need fine tuning of 
only potential parameters, the phantom field model also 
requires fine tuning of the initial condition of the field 
to account for the current acceleration of the universe. 
The model fits supernova la observations fairly well for 
a certain range of parameters. The best fitting param- 
eters of the model correspond to the equation of state 
parameter w<f, = —1.74 and can vary up to -2.4 at 95.4 
% confidence level. In our analysis we have investigated 
a particular model which avoids future singularity. It is 
expected that other models of this class (as mentioned 
above) would exhibit similar behavior. It would be inter- 
esting to further constrain the parameter space by other 
observations like CMB and structure formation. 
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